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Embedded Star Clusters in the W51 Giant Molecular 
Cloud. 
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ABSTRACT 

We present sub-arcsecond (0.35"-0.9"), near-infrared J,H,K band photometric ob- 
servations of six fields along the W51 Giant Molecular Cloud (W51 GMC). Our obser- 
vations reveal four new, embedded clusters and provide a new high-resolution (0.35") 
view of the W51IRS2 (G49. 5-0.4) region. The cluster associated with G48. 9-0.3 is 
found to be a double cluster enclosed in a nest of near-infrared nebulosity. We con- 
struct stellar surface density maps for four major clusters in the W51 GMC. These 
unveil the underlying hierarchical structure. Color-color and color- magnitude diagrams 
for each of these clusters show clear differences in the embedded stellar populations 
and indicate the relative ages of these clusters. In particular, the clusters associated 
with the HII regions G48.9-0.3 and G49.0-0.3 are found to have a high fraction of 
YSOs and are therefore considered the youngest of all the near-infrared clusters in the 
W51 GMC. The estimated masses of the individual clusters, when summed, yield a 
total stellar mass of ~10^Mq in the W51 GMC, implying a star formation efficiency 
of 5-10%. These results in comparision with the CO observations of the W51 GMC, 
suggest for the first time, that star formation in the W51 GMC is likely triggered by 
a galactic spiral density wave. 
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1 INTRODUCTION 

W51 is classically known to be a complex of compact ra- 
dio continuum sources representing a luminous star form- 
ing region in the galactic disk. The A21 cm maps show two 
major chunks of radio emission (representing complex HII 
regions) that are referred to as W51A and W51B. A dif- 
fuse and extended component of radio emission surround- 
ing and extending to the east of W51B is called W51C 
and identified as a supernova remnant. The HII region com- 
plexes W51A and W51B are associated with intense molec- 
ular emission that is spread over an area of l°xl°centered 
on (I,b)~(49.5°, -0.2°). This large scale molecular emission, 
which encompasses both W51A and W51B, is called the 
W51 G iant Molecular Cloud (GMO. ICarpenter fc Sanderj 
Jl998h Chereafter CS98) present CO J=l-0 and ^^CO J=l-0 
maps of the W51GMC and show that all of the bright, com- 
pact HII regions are situated along a high ve locity cloud 
which they identify as the "68km s"^" cloud. IKooI (1199911 
present higher resolution maps of the W51B region in CO 
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J=2-l and CO J=l-0 emission. CS98 argue that W51 is in 
the top 1% of all Galactic molecular clouds by size, and is 
in the top 5-10% by mass. 

To better describe this complex region in the con- 
text of the work presented here, we show in Fig.l a 
multi-wavelength view of the W51GMC. The top panel 
shows contou r s of t he A21 cm radio continuum image from 
IKoo fc MoonI 1119971) overlaid on a grey scale Sum emission 
image obtained by the MSX mission. The radio continuum 
sources (filled triangular symbols) G49.4-0.3 and G49.5-0.4 
are referred to as W51 A and the remaining components 
namely G49. 2-0.3, G49. 1-0.4 and G48. 9-0.3 constitute the 
W51B region. Mehringer (1994) used high resolution radio 
continuum and recombination line observations to identify 
several components of the HII regions associated with the 
W51 complex. Following his nomenclature, G49.5-0.4, the 
most luminous source in W51, is resolved into at least eight 
components (w51a-h). The most prominent among these are 
W51d and W51e which cor respond to the strong in frared 
sources IRS2 and IRSl dColdader fc Wvnn-WiUiamj|l99i) 
respectively. 

Fig. lb (bottom panel) displays ^^CO J=l-0 emission 
contours of the 68kms"^cloud (CS98) overlaid on a 2MASS 
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K-band image mosaiced by us. Comparison of Fig. la and lb 
shows that while the HII regions of W51A are fairly well en- 
closed by the observed ^ CO emission, the W51B HII regions 
extend beyond the densest regions of the molecular cloud, 
traced in 8//m emission. The mid-infrared MSX A-band im- 
age, which traces emission from the large-scale warm dust, 
resembles the 2MASS K-band image, which traces the em- 
bedded stellar content and hottest regions of gas and dust, 
in two ways. The extended, whitish regions in the 2MASS 
image (see also atlas image of W51 on 2MASS website), 
which are caused by extinction of the background light, co- 
incide with the brighter regions (darker shades) in the MSX 
image. However, these regions also coincide with compact 
clusters of stars in the 2MASS image. Further, these con- 
densations are also traced by the radio continuum emission, 
which delineates the HII regions and therefore the ionised 
gas. As is evident from Fig. lb, all of these bright conden- 
sations, which represent luminous star forming regions, are 
situated along the 68kms~^cloud. This indicates that they 
may be embedded in this giant molecular cloud. A cursory 
study of the 2MASS data showed that the regions enclosed 
by the bold, square boxes in Fig.l likely host embedded clus- 
ters. This prompted the observations discussed in this pa- 
per, where we present high spatial resolution near-infrared 
(NIR) photometric studies of the regions marked with the 
boxes in Fig. 1. However, this is not a complete sample, be- 
cause there are few other unstudied objects appearing on 
the 2MASS image that look quite similar to the regions of 
our choice. The following results provide a new picture of 
large scale star formation in the W51GMC, which was not 
clear from previous studies that mostly concentrated on the 
W51A region. 



2 OBSERVATIONS AND DATA REDUCTIONS 

Near-Infrared (NIR) observations were made at the 3.8 m 
United Kingdom Infrared Telescope (UKIRT) with the fa- 
cility imager UFTI. The [/FT/ pixel scale measures 0.091"; 
the available field of view is ~ 90" . Photometric observations 
through JHK broad band filters were secured for six fields 
during three nights of service observing (15 Sep 2000, 31 
May 2001 and 12 Aug 2001). Total exposure times in the J, 
H and K filters were 180, 180 and 360 seconds respectively 
on each of the target fields. The seeing varied between 0.32" 
to 0.9" with a mean of 0.5" in the ii'-band for the photomet- 
ric data. A nine point (3x3) jittered observing sequence was 
executed to obtain data that provided final mosaics with a 
total field of view of ^115" x 115". We note that the signal- 
to-noise ratio at the edges of these mosaics are lower than 
that within the central 90"area. Standard data reduction 
techniques involving dark subtraction and median-sky-flat- 
fielding was applied. 

We used tasks available under iraf for photometric 
analysis, daofind was used to identify sources in each im- 
age. A threshold equal to four times the mean noise of 
each image was found to be sufficient for this purpose. A 
psf model was computed by choosing several bright stars 
well-spaced out in the image. Photometry was performed us- 
ing DAOPHOT. Aperture corrections were determined by per- 
forming multi-aperture photometry on the psf stars. The in- 
strumental magnitudes were calibrated to the absolute scale 



Table 1 . Summary of observed fields 



Name 



Coords(J2000) 
RA DEC 



PSF Number of Stars 

"JHK JHK 



G48.9-0.3 19:22:14 14:03:09 0.36 662 1359 2071 593 

G49.0-0.3 19:22:26 14:06:46 0.50 954 1659 1757 702 

G49.2-0.3 19:23:02 14:16:41 0.90 469 882 905 293 

G49.4-0.3 19:23:05 14:28:07 1.10 542 702 744 304 

G49.5-0.4 19:23:40 14:31:07 0.32 651 977 1865 529 



using observations of UKIRT faint standards (FS29,FS35 
& FS149) that were obtained at airmasses closest to the 
target observations. The resulting photometric data are in 
the natural system of the Mauna Kea Consortium Filters 
jSimons fc Tokunagall2002h . For the purposes of plotting 
the se data in Fig.4 and F ig. 5 we have converted them to 
the ISessell fc Bretd lll988r (BB) system, since the main- 
sequence references are in the BB system. Since the direct 
transformation equations between the Mauna Kea system 
and the BB sytem are not established as of now, we first con- 
verted the Mauna Kea system to CI T system and then to th e 
BB system using equations given bv lHawarden et all i200ll) . 
We note that this indirect method can induce high errors for 
extremely red objects. Representative sub-images of G48.9- 
0.3 and G49.5-0.4, consisting of stars and nebulosity, were 
chosen for determining the completeness limits. Limits were 
established by manually adding and then detecting artificial 
stars of differing magnitudes. By determining the fraction of 
stars recovered in each magnitude bin, we have deduced 90% 
completeness limits of 18.2, 16.9 and 17.3 in the J, H and 
K bands respectively. Our observations are absolutely com- 
plete (100%) to the level of 16.2, 14.4 and 15.8 magnitudes 
in J, H and K respectively. The relatively higher detection 
level in K is due to the fact that the exposure times were 
twice that in the J k. H bands. The photometric accuracy 
of data presented here is generally <0.2mag although more 
than 30% of the data points presented in electronic tables 
have magnitude errors less than 0.05 mag. Astrometry for 
each of the images was calibrated by using the 2MASS as 
reference. The astrometric accuracy of the data presented in 
this paper is better than 0.5". 



3 EMBEDDED CLUSTERS IN W51 GMC 

We imaged six fields along the W51GMC as represented by 
the bold, square boxes in Fig. 1. The color plate in Fig. 2 
displays JHK band three color composites of each of the ob- 
served fields. Regions are referred to using the nomenclature 
of the associated UCHII region. The images are displayed on 
a linear scale except for the field G49.5-0.4 which is shown 
with a logarithmic scale to reveal the central region. As can 
be seen from Fig. 2, five of these six fields possess clusters 
of stars immersed in K-hanA nebulosity. The definition of 
a cluster here implies: (1) a significant increase in the stel- 
lar density with respect to an adjacent background, and (2) 
association with radio continuum sources. The field around 
G49.5-0.4 i s associated with the clusters know n as W51 IRSl 
and IRS2 JColdader fc Wvnn-WiUiamslll994D . The remain- 
ing four clusters are newly discovered from the observations 
presented here. The star symbols in Fig. 1 mark the positions 
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Figure 1. a) MSX A-band (8/im) image of tlie W51 GMC, displayed on a logarithmic scale, overlayed by 21 cm continuum contours 
JKoo fc Moonlll997r . The contour levels a re 0.1,0.2,0.4,0 .8 and 1.6 Jy beam~^. b) 2MASS A'-band image of the same region overlayed 
by contours of ^'^CO integrated intensity iCarpenter fc Sanders 1998^ . Contour levels are 8,16,32 and 64 K km s~^. The square boxes 
mark the observed regi ons, triangular symbo l s show the positions of associated UCHII regions and the star symbols represent the IRAS 
sources that satisfy the IWood fc ChurchwelJ ll989^ criteria. 
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Figure 2. JHK color composites of W51 Embedded clusters 



of IRAS sources that satisfy the lWood fc Churchwelil (Il989l) 
criteria for UCHII regions. Three of these IRAS sources are 
associated with embedded clusters; the remaining are unas- 
sociated with any NIR counterparts. The geometrical coin- 
cidence of the observed clusters with the ^^CO contours and 
the 21 cm contours indicate that these regions are "embed- 
ded" in the gas and dust of the W51 GMC. Table. 1 summa- 
rizes the positions and statistics for the five fields, derived 
from our JHK band photometric analysis. These numbers 



do not refiect the actual number of stars that belong to the 
clusters. Our photometric observations provide statistically 
meaningful data for only four of the five observed clusters; 
we therefore present detailed analysis only for these four 
clusters. The regions of G49.4-0.3 and G49. 1-0.4 are not in- 
cluded in the detailed analysis presented below. 
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3.1 Notes on individual clusters 

The four embedded clusters discussed below are associated 
with bright FIR and radio continuum emission. They are, in 
their own right, interesting objects for detailed high resolu- 
tion studies at various wavelengths. 

G49.5-0.4' This HIl region complex is composed of two 
major components known as W51d and W51e, which are as- 
sociated with two embedded clusters identified by the promi- 
nent infrared sources IRS2 and IRSl(see Fig. 1) The bright- 
est region, W51d, is well-studied at centimeter, millimeter 
and infrared wavelengths. W51d/IRS 2 is now known to be 
a young embedded star cluster that is larger and more lu- 
minous than the Orion- Trap ezium cluster . It contains the 
water maser complex W51N JForster fc Ca swell 19 891). sev- 
eral O H and SiO maser sources JCaswelJ i2003i' l iFish et alJ 
||2QQ3)) 8'iid is thought to contain several O type young 
stars. In Fig. 2, the inset of the G49.5-0.4 panel shows three- 
color JHK-composite image of the central region of the 
W51d/IRS2 cluster which, at a spatial resolution of 0.32", 
reveals more details than any other image of this region ob- 
tained until now. The region shows bright stars and neb- 
ulosity, the nebulous cloud appearing to be pinched off to 
form two lobes stretched east-west. These are called IRS2E 
and IRS2 W and are shown to h ave electron densities of 10^- 
10^ cm~^ l|Okumura et al.l2001rl . Our observations show that 
IRS2E is associated with a dense cluster, while IRS2W is rel- 
atively devoid of stars. Further, the high-resolution view sug- 
gests that IRS2E and IRS2W are probably two massive star 
forming clumps bound by a Roche-lobe potential. The cen- 
tral region is associated with 6-7 mid-infrared point-sources 
(;Kra eme r et al. 2001; Okamoto ct al. 2001), 6 of which ap- 
pear as point sources in our high-resolution NIR image. One 
of the sources appears as a red nebulous condensation (the 
brightest star of the wester n blob in the inset) a nd corre- 
sponds to the source KJD3 JKraemer et alJl200J) which is 
known to be associated with deep silicate absorption fea- 
tures. Most of these mid-infrared sources are thought to be 
O-type stars, although there is so me debate about whe ther 
they are 05-06 or 09 type stars JOkamoto et al.ll200lh . 

G48.9-0.3: This source is associated with one of the 
most spectacular clusters discovered in the W51GMC. The 
top right panel in Fig. 2 display the three color com- 
posite of this newly discovered embedded cluster. This 
cluster is located towards one edge of the ~100pc long 
68km s~^ molecular cloud. It appears to be placed on the tip 
of the molecular cloud and the observed HII region W51B. 
High resolution ISOCAM data for this region show strong 
FIR emission at 6/im with a morphology similar to the Sfim 
MSX data shown in Figl. This is a twin cluster immersed 
in a cradle of nebulosity visible at NIR wavelengths. The 
southern sub-cluster is a centrally symmetric cluster with a 
bright white/blue star situated at the geometric center. In 
contrast, the northern sub-cluster does not have a promi- 
nent source, although it does appear centrally symmetric. 
Indeed, here there appears to be a small ring of stars that 
is orientated in the plane of the sky. 

G49.2-0.3: This is the strongest UCHII region found in 
the W51B region and corresponds approximately to the lo- 
cation where the shock from the W51C supern ova remnant is 
thoug ht to be interacting with the W51 CMC JKoo fc MoonI 
Il9971) . It is associated with a group of bright red stars lo- 
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Figure 3. Stellar surface density contours generated by Nyquist 
binning K-band star counts. The dotted curves with grey shades 
represent the boundaries of the clusters and are obtained by us- 
ing a binwidth of 7". The levels of the dotted contours begins 
at 4 stars/pixel(7") and increase in steps of 2 stars/pixel. The 
solid curves represent the hidden structure of these clusters and 
are revealed by using smaller bin-widths represented by the grey 
circles in each panel. 



cated at the boundary of a cometary shaped HII region. 
The K-band (red) nebulosity has a similar morphology and 
is found to be well aligned with the radio free-free emission 
shown in Fig. 1. Unfortunately, our photometric data for this 
region suffers from relatively poor seeing (0.9") and a track- 
ing error which does not allow us to effectively resolve the 
otherwise dense group of faint sources that is enclosed by 
the cometary shaped HII region. 

G49.0-0.3: The bottom right panel of Fig. 2 shows 
the three color composite of this field. This region of the 
W51GMC is bright at 6-8/im and is relatively weak in terms 
of radio continuum emission, probably indicating a younger 
cluster and weaker HII region in comparison with other 
brighter sources in the W51B region. A faint "bow" of K- 
band nebulosity, that is similar in morphology to the associ- 
ated radio continuum emission, surrounds the cluster. This 
is probably a photo-dissociation region. Our three-color im- 
age shows a bright red star centrally-located in the cluster. 
Also found are two knots of K-band (red) nebulosity in the 
central region that can not be easily related to any particular 
star. 



3.2 Spatial distribution of stars 

It is evident from Figure, lb that the spatial distribution 
of near-infrared sources is non-uniform, on the global scale 
of the GMC and also on the local scale of the individ- 
ual clusters. In order to study the nature of this distribu- 
tion it is necessary to construct star count surface density 
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Figure 4. Colour-colour diagrams for four clusters. Tlie solid 
and broken lieavy curves represents the main-sequence dwarf and 
giant stars, respectively, and the dashed parallel lines are the red- 
dening vectors that enclose reddened main-sequence objects. The 
crosses on the reddening vectors are separate d by A^, =5 magThe 
dotted line indicates the locus of T-Tauri stars jMeyer et^L^^Q^) 
and t he dotted curve represents the HAeBe locus I Lada fc Adamsl 
Il992r . The star symbols in the bottom left panel represent the 
sources in the W51 IRS2 region. 



maps (SDMs). Such maps provide quantitative estimates 
of stellar surface densities and their variation across the 
clusters. They also reveal details about the nature of star 
formation and the hierarchical structure. Different meth- 
ods of generating such maps have been follo wed in various 
studie s of nearby young star clusters. While iLada fc Ladal 
il995h use overlapping rectilinear grids of boxes separated 
by Nyquist criteria to count the sources in the IC348 cluster, 
iHillenbrand fc HartmannI ((l99o) bin the source counts into 
10" pixels and smooth the resulting map with a symmet- 
ric gaussian profile with a 3-pixel width to produce stellar 
surface density maps of the Orion Nebula Cluster. 

In the above cases the targets of study were nearby 
star forming regions within 500 pc distance where meaning- 
ful source counts could be obtained over projected angular 
areas of several arcminutes. Massive star forming regions 
like W51 are generally at much large distances and thus 
subtend relatively small angular sizes making it difficult to 
study star count distribution. In these cases spatial resolu- 
tion plays an important role in the extraction of meaningful 
information from the SDM. The sizes of individual features 
we may expect to see on the SDMs depends on the Full- 
Width Half Maximum (FWHM) of the Point Spread Func- 
tion (PSF), the binning size adopted to generate the SDMs 
and the count level itself. Together, these determine whether 
or not a meaningful SDM can be made. The binning size is 
a smoothing factor that has to be carefully chosen in such 



a case. While a small bin size can include just a few sources 
and not provide a meaningful SDM, an overly large bin size 
can hide the true features in the SDMs. Given the high qual- 
ity of the seeing (0.32") in the A'-band images of two fields, 
we experimented with different bin sizes corresponding to 
6, 9, 12, 15, 18 and 22 times the FWHM of the PSFs. We 
found that while the majority of identifiable features did 
not vary significantly between 6 and 12 times FWHM, the 
features began to get smeared for a bin size larger than 12 
times FWHM. However, to identify the boundaries of the 
overall cluster and estimate its physical characteristics, a 
much larger bin size, to smooth all the features, had to be 
employed. 

In Figure. 3 we show SDMs for four fields where the de- 
tected source counts were sufficie nt to produce s uch a map. 
We have adopted the method of iLada fc Ladal (1 19951) who 
use a rectilinear grid of overlapping squares separated by 
half the width of individual squares to provide Nyquist sam- 
pling of the sources. We counted all sources from the A'-band 
photometry within our completeness limits discussed in the 
previous section. The cutoff K magnitude for this choice 
was 17.8 (85% detection limit); the magnitude errors were 
less than 0.3 mag. The dotted contours and grey shades in 
Fig. 2 indicate source distribution as obtained by a large 
smoothing box (150 pixels '^ 14"). These contours are used 
to identify the cluster boundaries and estimate the average 
surface density of each cluster. The thick contours overlayed 
on the grey shades represent SDMs generated by using a 
smaller binsize (14 x PSF) that shows smaller groups within 
these clusters. All the four clusters reveal a clumpy distri- 
bution of the NIR sources and the clumps are found to be 
either elliptical or circular in shape. The elliptical clumps are 
particularly noticeable in the maps of G48.9-0.3 and G49.5- 
0.4. The projected dimensions of these clumps are approxi- 
mately 0.1-0.2 pc. Using the SDM generated with large bins 
(smooth maps), we estimated the number of stars and the 
stellar surface densities for each cluster listed in Table. 2. 



4 COLOR-COLOR(CC) AND 

COLOR-MAGNITUDE(CM) DIAGRAMS 

Table. 1 lists the number of stars detected in each of the 
J,H and K photometric bands, and simultaneously in all 
three bands, for each of the fields. For the purposes of pho- 
tometric analysis of the clusters, we chose all those stars 
with photometric errors less then 0.2 mag. Tables. 2-5 list 
the astrometry(RA, DEC) and photometry (magnitudes and 
errors for each of the J, H and K bands) of the stars 
in sources G48.9-0.3,G49.0-0.3,G49.2-0.3 and G49. 5-0.4 re- 
spectively. For the four clusters with SDMs in Fig. 3., the 
stars detected in all three photometric bands are plotted 
on JHK color-color diagrams in Fig. 4. The solid and bro- 
ken heavy curves represents the main-sequence dwarf and 
giant stars, respectively, and the dashed parallel lines are 
the red dening vectors t h at en close reddened main-sequence 
obiects FOkumura et alj i2000h estimated the slope of red- 
denning vector in the direction of W51 to be ~1.9. Similar 
results are obtained by using 2MASS data. We therefore as- 
sumed a reddenning vector slope E(J-H)/E(H-K) = 1.9 (BB 
system) towards W51 region, which r elates to an interstellar 
reddening law with a value of R=3.12 iWhittet fc van Bredal 
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Figure 5. Color-magnitude diagrams for four fields. The vertical 
solid lines (from left to right) represent the main-sequence curve 
reddened by 0, 20 and 40 magnitudes respectively. The horizontal 
parallel lines represent reddening vectors. 

ll98(J). The dotted line indicates the locus of T-Tauri stars 
JMever et al.l |1997|) and the dotte d curve represents the 
HAeBe locus JLada fc Adanislll992ri . Stars that lie outside 
the region of reddened main-sequence objects are young stel- 
lar objects (YSOs) with intrinsic color excesses. By dered- 
dening the stars (on the CC diagrams) that fall within the 
reddening vectors encompassing the main sequence stars and 
giants, we found the visual extinction to each star. We de- 
redenned the stars to the K6-M6 part of the sequence of 
stars. The individual extinction values range from to 14 
magnitudes. From a histogram of these values, we estimate 
the average foreground extinction to each of the clusters. 
Two of the clusters (G48.9-0.3 & G49. 0-0.3) show two promi- 
nent peaks of extinction on the histogram corresponding to 
the two groups of points that can be identified on the CC 
diagram. The estimated values of extinction are listed in 
Table. 2. 

It can be noted from Fig. 4 that there are significant dif- 
ferences between the stellar population of the four clusters. 
G49.0-0.3 and G48.9-0.3 shows two distinct groups of points 
on the CC diagram separated by about 5-10 magnitudes of 
visual extinction. In contrast G49.2-0.3 and G49.5-0.4 do 
not show such distinct groups. Nearly 90% of the points 
that belong to G49.0-0.3 lie in the region of intrinsically 
reddened objects, thus signifying a high fraction of YSOs. 
Further, it appears that G49.0-0.3 has two groups of YSOs, 
one group more reddened than the other. We found that 
these two groups of YSOs are distributed homogeneously 
over the cluster area. Thus the more reddened group is likely 
embedded in the molecular clump while the least reddened 
group is probably partially out of or in front of the clump. 
The bright red star in G49.0-0.3 (see Fig. 2) at the center 



of the cluster shows colors similar to a T-Tauri star and is 
T^ highly reddened ( A„ ~35 mag) . This could be a massive star 
embedded in this cluster-forming core. 



In G48.9-0.3 there is a good mixture of both main- 
sequence and YSO populations. However, unlike G49.0-0.3, 
the two groups in this case are found to be distributed dif- 
ferently. The reddened group of stars are all found to oc- 
cupy the central parts of the two (northern and southern) 
sub-clusters (see Fig. 2), while stars without much redden- 
ing are found to be distributed over the entire FOV of the 
image. Therefore the average extinction value derived from 
the reddened group of stars is a good representation of the 
visual extinction to the molecular clump that hosts the two 
clusters. 



The G49.5-0.4 colour-colour diagram shows mostly 
main-sequence objects. There is a small percentage of YSOs 
with varying degrees of reddening. The star symbols in Fig. 4 
mark the positions of the stars associated with IRS2 (which 
is shown inset in the bottom-left panel in Fig. 2). These 
stars are thought to be 0-type stars. The placement of 
these sources in the upper-right part of the CC diagram 
supports the conjenctu re that these are indeed 0-type stars 
jLada fc Adam!Jll992^ . Finally, in Fig. 4 G49. 2-0.3 displays 
very few YSOs and only a few reddened main-sequence ob- 
jects. If we consider the ratio of sources that fall in the 
main-sequence zone and the YSO zone as an indicator of 
the cluster age, then the differences in the CC diagrams dis- 
cussed above imply that G49.0-0.3 is the youngest cluster 
in W51, while G48.9-0.3 is the second youngest cluster. In 
comparison, G49.5-0.4 and G49.2-0.3 are relatively evolved 
objects. 



Fig. 5 shows H-K vs K color-magnitude (CM) diagrams 
for all the four clusters. The vertical solid lines (from left to 
right) represent the main-sequence curve reddened by 0, 20 
and 40 magnitudes respectively. We have assumed a distance 
of 6.5 kpc for all the clusters, and respective A„ values (listed 
in Table. 1) to each cluster in order to appropriately repro- 
duce the main sequence data on this plot. Similar to the CC 
diagrams, the CM diagrams show clear differences between 
the four clusters. The reddened MS stars are found to be 
spread parallel to the MS line. Clusters G49. 5-0.4, G48.9- 
0.3 and G49.2-0.3 clearly displays this concentration of MS 
points, which is almost completely absent in the G49.0-0.3 
cluster. This agrees well with the results from the CC di- 
agram for G49.0-0.3, i.e. that YSOs comprise about 90% 
of the total membership of this cluster. The CM diagram 
of G48.9-0.3 shows the YSOs and MS stars well separated. 
Points that lie in the upper right corner of the CM dia- 
gram (above the line corresponding to 05) represent good 
candidates for massive YSOs. The relative number of such 
points are highest in G49.5-0.4, followed by G49.2-0.3 and 
G48.9-0.3. In contrast, G49.0-0.3 has almost all of its points 
concentrated in the lower right corner of the CM diagram, 
corresponding to a population of essentially low mass stars. 
However, an estimation of the true spectral type of each star 
can not be made using the CM diagram alone, and warrants 
NIR spectroscopic observations of the individual sources. 
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Table 2. Embedded Clusters in W51 CMC 





Number 


Stellar density (pc ^) 


Radius (pc) 


Age'^ 


Ai, 


Mass* 


Name 


of Stars 


Average 


Peak 


(\ 


/A/^) 


Myr 


mag 


Mq 


G48.9-0.3 


658 


175 


475 




3.43 


1.4 


(5,10) 


10000 


G49.2-0.3 


505 


130 


250 




3.48 


3.0 


(6,11) 


8000 


G49.4-0.3 


148 


123 


175 




1.95 


2-3 


5 


2000 


G49.5-0.4(IRS1D) 


309 


221 


535 




2.09 


0.7 


5 


4400 


G49.5-0.4{IRS1E) 


333 


179 


340 




2.39 


0.7 


5 


4700 



"^ age of associated HII region 

Calculated for an age of IMyr and A« ~7mag 



5 DISCUSSION 

5.1 Estimating Cluster Properties 

In Table. 2 we summarise some of the estimated properties 
of the four embedded clusters in the W51GMC. These esti- 
mates are based on the photometric analysis presented in 
previous sections and some standard methods. The clus- 
ters and their boundaries were identified using the SDMs 
shown in Fig. 3. The stars that fall within the outermost 
contours shown in Fig. 3 (dotted contours and grey shades) 
were counted to estimate the number of stars in the clus- 
ter. However, the stellar surface density (average and peak 
values) were estimated by using the star counts within the 
contour equal to half the peak value. It is clear from Fig. 3 
that the clusters are not circularly symmetric and indeed 
have irregular clumps. The projected area of the clusters 
were us ed to define an equivalent radius for the cluster 
as Yl=yJ A/-K, where A is the estimated area of each clus- 
ter from Fig. 3. We have estimated a mass for each clus - 
ter based on the method described bv lLada fc Ladal i2003l) . 
Following this method, we have assumed the K-luminosity 
function (KLF) mo dels of the Orion Trapezium Cluster 
JMuench et a l. '2002^ as a template and used the evolution- 
ary tracks of |DAntona fc MazzitcUi ( 1994) in computing the 
estimated masses. We have used an age of IMyr, distance of 
6.5Kpc and an average foreground extinction of A„~7mag 
to compute the masses. The effect of age on the mass of 
the cluster is such that ages ^l-2Myr yield lower masses 
compared to ages ^l-2Myr, with all other parameters fixed. 
However, these assumptions are fair, given the uncertainties 
involved in our assumed parameters. The estimated mass for 
the G49. 5-0.3 region (IRSl and IRS2) toge ther is QlOOMp) 
whic h is similar to the mass estimated by lOkumura et al.l 
i200Q) for this region (8200Mq) who assumes a Salpeter 
mass function. Our higher estimate is as expected because 
the spatial sampling here is a factor of 5-6 better compared 
tolOkumura ot al. (200(]j), thus revealing a larger number of 
cluster members. We also note that the estimated masses are 
quite sensitive to the assumed extinction values, particularly 
for distant sources such as these. Therefore, the estimated 
masses can be uncertain by a few thousand solar masses. 

The age of a cluster is one of the most difficult pa- 
rameters to estimate with our present knowledge because of 
bot h observational difficu lties and theoretical uncertainties 
fSee lLada fc Ladal ll2003t) for a discussion of this problem). 
In the case of massive star forming regions, some authors 



use the age of the associated HII regions as an indication 
of the age of the asso ciated groups of stars. For example, 
lOkumura et alJ (120001) give such ages for the four regions 
falling in the region of W51A. Extending this method to the 
entire W51GMC region, one finds that the G48. 9-0.3 is older 
than G49.5-0.4(IRS2). From the studies presented here, if we 
consider the ratio of YSOs to MS stars in each cluster as an 
indicator of its youth, and also consider the ages indicated 
by the CM diagrams, it turns out that the IRS2 region is 
older than the G48.9-0.3 region. This contradicts the result 
obtained by comparing the ages of the HII region. It has 
been suggested that the UCHII region phases of massive 
stars last for a time period much longer than lO ^yrs if they 
are born in dense and hot regions \de Preell995l) . This is be- 
cause the pressure exerted by the surrounding environment 
holds the expansion of the UCHII region phase. It may even 
be possible that massive star formation took place at a muc h 
later stage tha n the low mass star formation JHerbid il962l) : 
IStahleil (|l98a)). Therefore, neither the ages of the massive 
YSOs, nor the ages of the associated HII regions, are good 
indicators of the age of the associated embedded clusters. 



5.2 Hierarchical and triggered star formation in 
W51GMC 

The studies presented here indicate the existence of at least 
two, and possibly three, levels of hierarchy in the W51GMC. 
From Fig. 1 it is evident that star formation has occurred 
in two major clumps, known as W51A and W51B, within a 
~100pc-long molecular cloud. These clumps need not be real 
and could indeed be the effect of disruption by the super- 
nova remnant W51C (see lKoo fc Moor] (|l993) for a detailed 
discussion). The infrared image of the W51A region shows 
several clumps that are found as associations o f stars and 
nebulosity (see Fig. 1 of iHodapp fc Davia (|2002)). Next, as 
seen in Fig. 2, there is clustering within each of the four 
clusters. In particular, if we consider the example of G48.9- 
0.3, which is a double cluster, there appears to be an addi- 
tional level of hierarchy. All of the clusters, at various lev- 
els, indicate hierarchical clustering. Hierarchical structure 
is believed to be a direct consequence of turbulent energy 
th at operates simila r ly at diff erent scales (See, e.g. , reviews 
bv lElmegreen et al.l J2000l) fc iMac Low fc KlessenI J2004l) ). 
Thus the observed hierarchical clustering indicates the pres- 
ence of significant turbulent energy in W51GMC which is 
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not yet dissipated. The observed star formation has there- 
fore taken place in the presence of strong turbulent energy. 

If we sum-up the mass of all of the clusters listed in 
Table. 2 we find a total stellar mass of ~3O,OOOM0 . Our ob- 
servations probably represents only ~50% of the total star 
formation in the W51GMC, due to various reasons. Firstly, 
we have not sampled all the smaller associations/clusters 
that can be seen in a large scale 2MASS K-band image. 
Next, our observations have poor spatial resolution for two 
of the major clusters; this results in an overly undersampled 
cluster membership for each cluster. Also, the distributed 
component of individual young stars across the cloud can 
contribute a significant percentage to the stellar mass which 
is neglected. Thus, the total mass estimates from these 
clusters alone are lower limits. In view of these uncertain- 
ties, it may not be unreasonable to assume that the over- 
all star formation in the W51 GMC is associated with a 
mass of ~5-6x1O''M0. In such a case, the estimated mass of 
the W51GMC, -lO^M©, derived from observations of the 
molecular cloud itself (CS98), implies a star formation effi- 
ciency (ratio of stellar mass to cloud mass) of ~5-6%. CS98 
showed that four of the radio continuum sources, namely 
G48.9-0.3, G49.0-0.3, G49. 1-0.4 and G49. 2-0.3, are spatially 
and kinematically associated with the 68kms~^cloud alone. 
If we believe that all the observed clusters in W51 GMC 
do indeed belong to the 68kms~^cloud, then the mass es- 
timate of ~2x10'^Mq for the 68kms"^cloud (CS98) , to- 
gether with the mass estimates for the clusters in Table. 2, 
indicates star formation efficiencies ~10-15%. This value 
is comparable to the star formation efficiencies in some 
relatively nearby massive star forming clusters, such as 
NGC6334 (25%)(Tap i a et aL.199&) and W3 IRS5 (6-18%) 
jMeeeath et alJll996D . 

Low mass star formation in the solar neighbourhood is 
known to have a low efficiency, of the order of 1-3%. Trig- 
gering mechanisms are thought to result in the formation of 
clusters; they are also thought to be associated with higher 
star formation efficiencies. In the W51GMC, the embed- 
ded clusters are distributed inside a ~100pc-long cloud. The 
clusters also likely have ages of 0.5-3 Myrs because that is 
the commonly observed ages of observed embedded clusters 
JLada & Lada 2003). This suggests near-simultaneous star 
formation across the entire cloud, an event which could have 
been spontaneous and/or triggered by an external mecha- 
nism. If a triggering mechanism was involved then it would 
have to have been a large scale shock wave to influence the 
100 pc-long region. Internal trigger mechanisms are very un- 
likely in this case, since the crossing time for the cloud, for an 
average velocity of ~10kms~^, is ~10Myr, which is larger 
than the typical life time of 3-5 Myr for the embedded clus- 
ters. When we consider these facts and note the location of 
the W51GMC at t he tangent point of the Sagittarius spiral 
arm iBurtonlll97(J) . it seems likely that the observed star 
formation in the W51GMC is triggered by galactic spiral 
density waves. Giant clouds like W51 are indeed thought to 
undergo gravit ational collapse u nder the influence of spiral 
density waves llElmegreenlll994r) . Therefore, W51GMC may 
be one of the best examples in our galaxy of triggered star 
formation by spiral waves. 



6 SUMMARY AND CONCLUSIONS 

We present sub-arcsecond resolution, near-infrared J,H,K 
band photometric observations of six flelds in the W51 GMC 
(which includes the W51A and W51B regions). We con- 
struct a multi- wavelength view of the W51 GMC, combining 
CO and A21cm observations from the literature with MSX 
and 2MASS data. All of these observations together suggest 
a unified view for W51, the HII regions being associated 
with a ~100pc long molecular cloud named the "68km s~^" 
cloud by C arpenter fc Sanders ( 1998) . The results of our 
near-infrared observations can be summarised as follows. 

l)From our observations of six regions, four new, em- 
bedded clusters are discovered that are associated with the 
UCHII regions G48.9-0.3, G49.0-0.3, G49.2-0.3 and G49.4- 
0.3. Among these, the cluster associated with G48.9-0.3 is 
found to be a twin cluster enclosed in a nest of near-infrared 
nebulosity. Our observations also provide a sub-half arcsec- 
ond view of the W51 IRS2 (G49.5-0.4) region that is known 
to be forming several O-type stars. This high resolution im- 
age provides a view which suggests that the W51 IRS2E and 
IRS2W are two massive star forming clumps that are likely 
connected by a Roche-lobe potential. 

2) Using K-band star counts, we construct stellar sur- 
face density maps for four clusters, (G48.9-0.3, G49.0-0.3, 
G49.2-0.3 and G49.5-0.4.) These maps unvefl the underlying 
hierarchical structure and allow us to estimate some physical 
parameters for the clusters. Color-color and color-magnitude 
diagrams for each of the clusters indicate distinct differences 
in their embedded stellar population. The clusters associ- 
ated with G48.9-0.3 and G49.0-0.3 are found to contain a 
high percentage of YSOs, suggesting their relative youth in 
comparision to the other regions studies. 

3) Assuming an Orion Trapezium IMF and theoreti- 
cal pre-main sequence tracks, we estimate the masses of the 
individual clusters. The combined mass of all the clusters 
amounts to a net stellar mass of ^lO'^M© in the W51 GMC. 
This estimate, when combined with est imates of the molec- 
ular cloud mass from lCarpenter fc Sanders (_1998j), suggests 
a star formation efficiency of ~10%. 

Finally, the combined view of the W51 GMC, from our 
new observations and from the literature, suggests that star 
formation in the W51 GMC is probably triggered by a galac- 
tic spiral density wave. Thus, W51 may represent a relatively 
unique example of massive star formation in our galaxy. 
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